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The resistivity minima of the quantum HaU effect arise due to the locahzation of the electron 
states at the Fermi energy, when it is positioned between adjacent Landau levels. In this paper we 
calculate the localization length ^ of such states at even filling factors v = 2N. The calculation is 
done for several models of disorder ("white-noise," short-range, and long-range random potentials). 
We find that the localization length has a power-law dependence on the Landau level index, ^ cx A'^" 
with the exponent a between one and depending on the model. In particular, for a "white- 
noise" random potential ^ roughly coincides with the classical cyclotron radius. Our results are in 
reasonable agreement with experimental data on low and moderate mobility samples. 



I. BACKGROUND AND RESULTS 

The appearance of narrow resistivity peaks separated 
by deep minima is a defining feature of the quantum Hall 
effect (QHE).EJ The explanation of such a dependence of 
Pxx on the magnetic field B is based on the idea of local- 
ization. The states at the Fermi energy are localized at 
almost all B except for a few discrete values Bat where 
the Fermi energy is at the center of A^th Landau level 
(LL). Near such special values the localization length ^ 
is believed to diverge. 
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where 7 is a critical exponent. The analytical calcula- 
tion of 7 is a notoriously difhcult problem. (Numeri- 
cal methods give 7 = 2.35 ± 0.03, see Ref. I). At the 
same time, the calculation of ^ away from the critical re- 
gion turns out to be much simpler. Such a calculation 
is the subject of the present paper. As a demonstration 
of the method, we calculate ^ at discrete values of B, 
B ~ {Bn + i?Ar_|_i)/2. They correspond to the minima of 
Pxx in the transport measurements. 

Generally speaking, the definition of ^ is not unique. 
In this paper we will adopt the following one. 
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where ip{r) is the wave function of the state at the Fermi 
level. The averaging is assumed to be done over the dis- 
order realizations. 

Our definition of the localization length is chosen to 
represent an experimentally measurable quantity. In- 
deed, it is well known that transport at sufficiently low 
temperatures proceeds via the variable-range hopping. 
In its turn, the hopping conduction is determined by the 
typical decay rate of the tails of the wave functions. Def- 
inition (1.2) relates ^ precisely to this typical rate. 

Let us further elaborate on this point. The bulk of 
low-temperature experimental data on the quantum Hall 



devicesi~il can be successfully fit to the following depen- 
dence of Pxx on the temperature T, 
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which can be interpretecB in terms of the variable-range 
hopping in the presence of the Coulomb gap.Q In this 
theory Tq is directly related to ^ defined by Eq. (1.2), 
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where e is the electron charge and k is the dielectric con- 
stant of the medium. Using Eq. ( |1.4| ), one can extract 
the dependence of ^ on i? from the low-T transport mea- 
surements in the straightforward way. (We will discuss 
experimentally relevant issues in more detail in Sec. 0.) 

In this paper we calculate ^ using a model where the 
disorder is described by a Gaussian random potential 
U (r) with the two-point correlator 
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We will assume that function C(r) becomes small at r 
larger than some distance d and that C (r) does not have 
any other characteristic lengths. The rms amplitude of 
the potential, ^/C(0), will be denoted by W. We will as- 
sume that W is much smaller than Ep , the Fermi energy. 
The electron-electron interaction is ignored at this stage. 

To facilitate the presentation of our results we would 
like to introduce the phase diagram shown in Fig. |l|. The 
vertical axis stands for the dimensionless parameter kpd, 
where kp is the Fermi wave- vector of the two-dimensional 
electron gas, kp = \phrn, n being the electron gas den- 
sity. The horizontal axis is the LL index N ~ (^^0^/2, 
where I = y^hc/eB is the magnetic length. The Fermi 
level is assumed to be at the midpoint between the cen- 
ters of iVth and TV -I- 1st LLs. The axes are in the loga- 
rithmic scale. The ratio W/Ep is assumed to be fixed. 
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FIG. 1 

monoscale random potential [Eqs 

entire phase space is divided into regions with different 
dependence of ^ on A'^ and kp d. T he boot-shaped region 
AHIGEB is described by Eq. ( 1.13 ), the shaded region to 
the right of the hne DFGI by Eqs. (0) and the re- 

gion BEGFD by Eq . ([l.iol ), and finally, the region below 
the line HI by Eq. ( |l.l5| ). The arrows show the "trajec- 
tory," traced by a "standard" sample (see definition in 
Sec. |v| ) as the magnetic field decreases. 

Several lines drawn in Fig. ^ divide the phase plane 
into the regions with different dependence of ^ on _B. Let 
us explain the physical meaning of these lines. The line 
BEGI is the line where the densities of states of neigh- 
boring LLs start to overlap as N increases {B decreases) . 
Thus, to the right of this line the density of states at the 
Fermi level is practically equal to its zero field value. We 
will refer to this region of the parameter space as the re- 
gion of overlapping LLs. To the left of the line BEGI only 
the tails of the neighboring LLs reach the Fermi level and 
the density of states is much smaller than at zero field. 
This region will be called the region of discrete LLs. The 
equation, of the line EGI is (cf. Refs. ^|-|lo|) 
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Another line in Fig. GFD, separates the regions of 
different dynamic properties. To the left of this line 
the guiding centers of the cyclotron orbits would per- 
form the regular drift along certain closed contours. This 
phenomenon has been dubbed "classical localization" in 
To the right of the line GFD (shaded sector in 
the motion of the guiding center is diffusive (on 
not too large length scales). The equation for the line 
GFD has been derived in Ref. [ill. 
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In the diffusive region the calculation of ^ reduces to 
the calculation of the "classical" conductivity a^x by 
means of the ansatz (for discussion and bibliograthy see 
Refs. I and ^ 
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The "classical" a^x is to be calculated by virtue of the 
Einstein relation, i.e., as a product of quantum density of 
states and the classical diffusion coefficient. The phys- 
ical mechanism of the localization in this region is the 
destructive interference of the classical diffusion paths. 
The calculation of the "classical" axx to the right of the 
line GFD and in the logarithmically narrow sector to its 
left (where dxx is still larger than ^) has been do ne in 
Ref. Oin some detail. As one can see from Eqs. ( |l.6| ) 
and ([l77| ), the studied region corresponds to rather long 
range of t he r andom potential, kpd > {Ep /W)"^^^ . How- 
. Eq. (P 



ever, 
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I applies for smaller values of parameter kpd 
as well. That is as long as we stay to the right of the line 
GI. In the entire shaded sector to the right of DFGI the 
motion is diffusive. The corresponding "classical" axx is 
given by the usual Drude-Lorentz formula 
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Equations ( |1.8| ) and ( |1.9D enable one to calculate ^ up to 
a pre-exponential factor. In this sector ^ is exponentially 
large. 

As N decreases and the boundary DFGI of the diffu- 
sive region is crossed, the "classical" axx rapidly falls off. 
Above the point G this is brought about by the afore- 
mentioned "classical localization"; below the point G it 
is caused by the rapid decrease of the density of states 
at the Fermi level. Already slightly to the left of the line 
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DFGI the " clas sical" axx becomes much less than ^ and 
the ansatz (1.8) loses its domain of applicability. On the 
physical level, the nature of the particle motion changes: 
the diffusion is replaced by quantum tunneling. This pa- 
per is devoted to the calculation of ^ in the tunneling 
regime (unshaded area in Fig. |l|). 

One has to discriminate between the tunneling in 
the case of overlapping LLs and in the case of discrete 
LLs. The former is realized in the region above the line 
BEGFD. The idea of the derivation of ^ in this regime 
belongs to Mirnikov_and Sokolovll3 who applied it to the 
lowest Landau level£3 The argument goes as follows. In 
the described regime the density of states near the Fermi 
level is high. On the quasiclassical level such states can 
be thought of as a collection of close yet disconnected 
equipotential contours, along which the particle can drift 
according to the classical equations of motion. Nonzero 
^ appears as a result of the quantum tunneling through 
the classically forbidden areas between adjacent contours. 
The localization length is determined by the spatial ex- 
tent of relevant equipotential contours and by the char- 
acteristic tunneling amplitude. 



Reference 12 has been criticized in literature for, e.g., 
neglecting the interference effects. In our opinion this 
criticism is unjustified. The authors of Ref. na have 
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clearly indicated the domain of applicability of their the- 
ory. It can be verified that within this domain the ampli- 
tude of tunneling between the pairs of contours is small; 
hence, the probability of return to the initial point after 
at least one tunneling event is also small. In this case the 
interference phenomena can be safely ignored (cf. Ref. ||). 

The case of high Landau levels requires some modifi- 
cations to the original method of Mil'nikov and Sokolov. 
The details are given in Sec. 0. We have found that the 
region of the phase space bounded by the line BEGFD 
consists, in fact, of three smaller regions (see Fig. |l|) with 
different dependence of ^ on A^: 



valid in a much larger domain. The differences between 
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(the equation labels match the region labels in Fig. |^). 
The boundary EC, which separates the regions I and II, 
is given by 



EC : kpd - N. 



(1.11) 



The variety of different subregimes in the region BEGFD 
appears because of an interplay among three important 
length scales of the problem: the correlation length d 
of the random potential, the cyclotron radius Rc = 
{2N + l)/kp, and the percolation length ^pcrc (the typical 
diameter of the relevant equipotential contours) . In this 
connection note that Eq. ( [l.ll ) is simply Rc ~ d. 

Compared to such a variety, the situation in the region 
of discrete LLs (AHIGEB) is very simple: the depen- 
dence of ^ on N is given by a single formula. Suppose 
that the Fourier transform C{q) of the correlator C(r) 
[see Eq.(O)] has the form 



C{q) = C(0) exp 
then ^ is given by 
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(1.13) 



^^-1 W J d 
The logarithmic factor neglected, this can be written as 
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Formula (1.13) has previously ajapeared (for /3 = 2) in the 
work of Raikh and ShahbazyanJlil These authors consid- 
ered the case of the lowest LL {N = 0) but suggested that 
it is also valid for Rc <C d, i.e., within t he kn ife-shaped 
region AHEB. We demonstrate that Eq. (1.13) is in fact 
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this work and Ref. |lj are outlined in the end of Sec 

The only part of the phase diagram we have not dis- 
cussed yet is the area of "white-noise" potential. It is 
located below the line HI, i.e., kpd = 1. The correspond- 
ing formula for ^ is 
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which match es Eq. (1.12) at kpd ^ 1. Previously, fcp^. 
mula ( 1.15b ) has been obtained by Shklovskii and Efrostj 
and also by Li and ThoulessE3 for the lowest LL, = 0. 

Ne glecti ng the logarithmic factor, we can write 
Eq. (1.15a) in a remarkably simple form. 



(1.16) 



Remarkably, the quantum localization length ^ is deter- 
mined by a purely classical quantity, the cyclotron radius! 

Th e ba sic idea used in the derivation of Eqs. (1.12) 
and ( 1.15 ) is to study not the tunneling of the particle 
itself but the tunneling of the guiding center p of its cy- 
clotron orbit. For definiteness, consider the tunneling in 
the y-direction. The effect of the magnetic field can be 
modelled by means of the effective "magnetic" potential. 



(1.17) 



acting on the particle. The classical turning points for 
this type of potential are at the distance Rc from py. 
Therefore, if py does not change its position, the longest 
distance that the particle can travel without getting un- 
der the magnetic barrier is 2Rc. And since the barrier 
increases with y, the suppression of the wave function, 
which starts beyond this distance, is faster than a simple 
exponential. 

In the absence of the random potential, py is a good 
quantum number; however, if the external potential is 
present, it can scatter the particle, which would cause a 
change in the guiding center position. Such a "scattering- 
assisted'L-tpmcling modifies the overall decay of the wave 
functionBEllla from the super-exponential to the plain 
exponential one. 



^-(0,2/) 



Denote a typical displacement of the guiding center af- 
ter one scattering act by Apj,. The physical picture of 
tunneling depends on the relation between Ap^ and 2Rc. 

The case Apy > 2i?c, which is typically realized at the 
lowest LL, has been studied previously in Refs. 15- 1^. 
In this case the tunneling involves the propagation un- 
der the magnetic barrier. Note that the barrier itself 
no longer increases as y squared, which would be with 
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Py = const [see Eq. (1.17)]. After a series of displace- 
ments of the guiding center, the barrier acquires a saw- 
tooth shape. In this regime the under-barrier suppression 
is an important factor in the overall decay of the wave 
function. 

In contrast to the lowest LL, at high LLs {N ^ 1) 
where Rc is large, the inequality of the opposite sense, 
i.e., Apy < 2Rc, is typically realized. In this case the 
particle does not propagate under the magnetic barrier 
at all! However, the tunneling a distance y ^ Apy 
requires a large number M y/Apy of the scatter- 
ing acts. The amplitude of each act is proportional to 
W and is inversely proportional to a large energy de- 
nominator E ^ Ep — hLUc{N + ^). At the resistivity 
minima of the QHE, which we are mainly interested in, 
E = TiuJc/'2 ^ W (discrete LLs), which implies that 
the typical scattering amplitude is small and that the 
wave function decays exponentially with y even though 
the electron never propagates under the magnetic bar- 
rier (this argument is simply a verbal representation of 
the locator expansion). 

The case of the "white noise" potential is quite illu- 
minating in this respect. The optimal tunneling path is 
sketched in Fig. |^. The optimization is based roughly on 
the requirement that each scattering event should dis- 
place the guiding center by the largest possible distance 
without placing the particle under the magnetic barrier. 
Clearly, this distance is equal to 2Rc, which makes the 
tunneling trajectory look like a classical skipping orbit 
near a hard wall, see Fig. |[ 
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FIG. 2. The optimal tunneling trajectory of an elec- 
tron in the random potential of the "white-noise" type. 
Crosses symbolize the scattering acts. The direction of 
the y-axis is from the left to the right. 



Let us estimate the localization length corresponding 
to this optimal path. Since Apy ~ 2i?c, the number 
of the scattering events needed to travel the distance y 
is M ~ y/{2Rc). As discussed above, after each event 
the wave function decreases by a factor of the order 
of W/{Tiu)c)- Hence, the overall suppression factor is 
(?iWc/VF)~*^, which means that ^ ^ 2Rc/ \n{fiuJc/W) in 
agreement with Eq. ( 1.15a ). This derivation will be done 
more carefully in Sec. III . 

Formula (1.13) can be derived in a similar way. Af- 



ter each scattering event the wave function decreases by 
a factor of the order of / = {W/ E)ex.^[- {1/2 13) [qdf] 
where q — Apy/P is the typical wave vector absorbed 
in the scattering act. The total suppression factor after 
propagating the distance y is of the order of / to the 
power yj Apy. Optimizing this supp ressio n factor with 
respect to Apy, one arrives at Eq. ( 1.13| ). A detailed 



derivation will be done in Sec. III. 

The paper is organized as follows. Section |l| is de- 
voted to general considerations and qualitative deriva- 



tion of Eq. ( |1.13D and (1.15). In Sec. |HI| this derivation is 
made more rigorous assuming that the random potential 
is short-range (or "white- noise" ) . In Sec. |^we consider 
the long-range potential in the regime of discrete LLs. 
The approach is diffe rent from that of Sec. |l^ but the 
final result. En. (1.12), is the same as for the short-range 
case. In Sec. ^ we consider t he ca se of overlapping LLs. 
The variety of regimes in Eq. (1.10) is explained with the 
help o^jresults developed in the field of statistical topog- 
raphy.ll3 Finally, in Sec. VI we compare our results with 



available experimental data for moderate mobility sam- 
ples and propose a method to perform the measurements 
with modern high-mobility devices. 



II. GENERAL CONSIDERATIONS 

Following the overwhelming majority of papers in the 
field, we will take advantage of the single Landau level ap- 
proximation. In this approximation the original Hilbert 
space is truncated to the functions, which belong to the 
particular (A^th) Landau level. It is conventional to 
choose the orthonornial set of functions 



V'„(a;,j/) = ^=$Ar(2/ + p„), (2.1) 



V2^n 
2TrP 



(2.3) 



-n, n = 0, 1, 



to be our basis states. Here Lx is the a;-dimension of 
the system and Hn(z) is the Hermite polynomial. Such 
functions are the eigenfunctions of the Hamiltonian 



(2.4) 



in the Landau gauge, A — {By, 0,0), provided there is 
no external potential {U = 0). 

The single Landau level approximation works well if 
the cyclotron frequency ujc is the fastest frequency in the 
problem. This is obviously the case for the discrete LLs, 
i.e., in the region AHIGEB in Fig. |l|. It is less trivial and 
it was demonstrated in Ref. ^ that the inter-LL transi- 
tions are suppressed in the region above the line BEGFD 
as well. When such transitions are neglected, the guiding 
center coordinates, py = y—{vx/oJc) and px = x+{vy/u!c) 
become the only dynamical variables in the problem. 

Since the random potential is assumed to be isotropic, 
so is the ensemble averaged decay of the wave functions. 
With the above choice of the basis, however, it is conve- 
nient to study such a decay in the y-direction: from the 
point (0,0) to the point (0,2/). 
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As the next step we notice that the guiding center co- 
ordinates satisfy the commutation relation 

[PyiPx] = il^- 

Thus, py plays the role of the canonical coordinate while 
the quantity {h/P)px is the canonical momentum. It 
is therefore natural to use the py-representation for the 
wave functions. For example, in this representation 



wave functions (2.1) become delta-functions. In general, 
the transformation rule between the two representations, 
ip{x,y) and (t>{py), is given by the formula. 



dPy 



'^N{y~Py)4>{Py)e 



(2.5) 



Definition (1.2) of the localization length can also be 
written in terms of the electron's Green's function G, 

I 



hm -(ln|G(0,t/;i?)|) 

y^oo y 



(2.6) 



On the other hand, Eq. ( |2.5| ) leads to the following 
Green's function transformation rule for G(0, y; i?), 



G= /^g^$jv(pi)<i>iv(p2)Gp(pi,p2 + 2/;i?), (2.7) 



where Gp(pi, p2 \ E) is the Green's function in the guiding 
center representation (as in Sec. ^ the energy E is refer- 
enced to the Landau level center, and so we are interested 
mostly in the case \E\ = fiujc/2). 

Green's function Gp{pi,p2',E) satisfies the 
Schroedinger equation (with the delta-function as a 
source) 



Gp{pi,p2;E) 



Hpi - P2) 
E 



-E 



Uo{qx,Pl + -qxl^)Gp{pi + qj'^,p2;E), 



(2.8) 



where the tilde indicates the Fourier transform over the 
corresponding argument. The quantity Uq has the mean- 
ing of the random potential averaged over the cyclotron 
orbit (cf. Refs. 



[compare with the definition of ^, Eq. (2.6)]. Unfortu- 



Uoiq) 



FNiy) = En [t^ 



U{q)FN{qn, 



-yVii' 



(2.9) 
(2.10) 



Ln{z) being the Laguerre polynomial (the tilde over the 
symbol itself indicates the two-dimensional Fourier trans- 
form). 

As discussed in the previous section, in the absence 
of the random potential, py is a good quantum num- 
ber. The tunneling requires propagation under the mag- 
netic barrier, which leads to a super-exponential decay of 
G{0,y;E). Indeed, in the absence of the random poten- 
tial Gp{pup2;E) = 5{qi_^ P2)/E [see Eq. Upon 



substitution into Eq. (2.7) one recovers the well-known 
expression for the Green's function in the clean case. 



G{0,y;E) 



FN{y) 

2nPE ' 
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Since i^Ar (y) i s the product of a polynomial and a Gaus- 
sian [see Eq. ( 2.10 )], the Green's function G{0,y;E) de- 
cays faster tha n the exponential at large y. In view of 
definition ([2.6^ this means that ^ = 0. In fact, the struc- 
ture of Eq. (^.71) suggests that nonzero ^, i.e., a simple 
exponential decay of G, is possible only if Gp decays no 
faster than a simple exponential. In other words, ^ is 
nonzero only if ^p is nonzero where 



1 = - lim -{\n\GpiO,p;E)\} 



(2.12) 



nately, ^ and may differ. Only the inequality 



(2.13) 




FIG. 3. Tunneling paths of the guiding center 
(schematically). Green's function G(0, y; E) is a sum over 
the paths near = 0; Green's function Gp{Q,y; E) is a 
sum over all the possible paths connecting the lines py — 
and Py = y. 



Indeed, Gp{pi,p2]E) typi- 

-\P2-pi\/ip+i4>(pl,P2) ^ Jf ^l^g 



is guaranteed to be met. 
cally behaves like Gp ~ 
phase (j){pi, P2) is a smooth function of coordinates, then 
|G(0,y;£;)| - \Gp{fl^]E)\, and C = ip- Otherwise, the 
integrand in Eq. ( |2.7| ) oscillates rapidly, \G{Q,y\E)\ <C 
|Gp(0, y; E)\ and ^ < ^p. On a physical level, ^ describes 
the tunneling between two point-like contacts while ^p 
characterizes the tunneling between two infinite paral- 
lel leads. While the Feynman paths contributing to the 
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former process make up a narrow bundle near Px = ^ 
(Fig. H), the latter one gathers contributions of many 
such bundles. As a result, the amplitude of the latter 
process is much larger on the account of rare places ( "pin- 
holes") where the tunneling is unusually strong. 

Nevertheless, in many cases f and are very close to 
each other; for instance, when the random potential is 
short-range (see the next section). 

Next we would like to present a simple yet very in- 
structive model. This model has a great advantage of 
being solvable. 

Suppose that the averaged random potential has the 
form 



C/o(p) = 



U;{Py), (2.14) 



where Ui, ReC/2, and Imf72 are mutually independent 
Gaussian random variables. Similar to the above, we 
will assume that they have amplitude W and correla- 
tion length d. A simpler model with pj^-independent Ui, 
KeU2, and Im[/2 has been studied in Rcfs. |l^ and |l9|. 

For the model potential (2.14) all the points along the 
Px axis are statistically equivalent; the "pinholes" we 
mentioned above are absent; therefore, ^ = £,p- 



we see 



From Eq. 
{pi\Uo\p2) is zero unless pi 



that the matrix element 
= P2 or pi~ p2 — ±qP. It is 
convenient to assume that qP is divisible into 2TrP/Lx, 
the smallest dist ance between the centers of gravity of 
the basis states (2J). In other words, we will assume 
that Lxq/{2TT) is an integer. Under this condition, the 
system can be split into Lxq/{2'K) independent chains 
(Fig. 1^). The guiding center coordinates {pn} in each 
chain form an equidistant set: Pn+i ~ Pn = qf. The 
hopping is allowed only between the nearest neighbors 
of the same chain and is characterized by the hopping 
amphtude C/2 + Pn+i)/^]. As for Ui{pn), it plays the 
role of the on-site energy. 







qP 
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FIG. 4. The model system. Each vertical tick corre- 
sponds to one of the basis states -i/in. The distance between 
the ticks is 2-k]} jL^- The hopping (symbolized by arrows) 
is possible at (much larger) distance gZ^, and takes place 
between the states forming an equidistant chain. Two of 
such chains (one marked by the dots and the other by the 
triangles) are shown. 

The localization length of a disordered chain is given 
by the exact formula due to Thouless,E3 which in our case 
reads 



V= \dE'D{E')\^\E-E'\-{\^\U2{p)\). (2.15) 



D{E) being the disorder- averaged density of states nor- 
malized by the condition JdE D(E) = 1. It is notewor- 
thy that D{E) can in principle be found exactly if all 



the matrix elements are statistically independently i.e., 
if qP' ^ d. If this is not the case, then D{E) can be 
calculated by some approximation scheme. At argi-rate, 
D{E) is small if \E\ > W (see e.g., Lifshitz et alH) In 
this case we can expand the logarithm in Eq. (2.15) in 
the powers of E' /E to obtain 
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(2.16) 



Taking the average in Eq. (2.16), we obtain 
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(2.17) 



where C = 0.577 ... is the Euler constant. 

It is quite remarkable that ^ does not depend on 
whether the successive hopping terms U2[{pn + Pn+i)/2] 
are correlated or not. With the high degree of accuracy, 
0{W'^ / E"^), the localization length has the same value 
for qP > d ("short-range" disorder) and qP < d ("long- 
range" disorder). 

The qualitative derivation of Eq. ( 2.17 ) can be done 
with the help of the locator expansion (see a similar ar- 
gument in the preceding section) . The tunneling through 
a distance y is achieved by a minimum of M = y/{ql'^) 
hops. Each hop is characterized by the hopping ampli- 
tude of the order of W/E. Thus, the suppression factor 
of the wave function over a distance y is of the order 
of iW/E)^ . On the other hand, this factor is equal to 



which leads to Eq. (2.17) 



Let us now return to the original problem with t he tw o- 
demensional random potential [Eqs. ( |l.5| ) and (1.12)]. 
Leaving a more rigorous calculation for the next two 
sectio ns, w e will pres ent heuristic arguments leading to 
Eqs. (Ml ) and ([lTsI ). 

Let us divide the entire spectrum of Fourier harmonics 
of the random potential Uq into bands q„ — ^Ag < qx < 
qn + jAg, n = 1, 2, ... of width Aq ^ 1/d. Denote by 
Uq{p) the combined amplitude of the harmonics, which 
make up the band centered at g. 



dqx 
27r 



Uoiqx,Py) 



(2.18) 



q~h^q 



If g 3> l/d, the corresponding band is very narrow, and 
Uq as a function of px looks very much like a plain wave, 
Uq oc 6*'^''== , exactly as U2 in the model problem. Sup- 
pose that the scattering caused by different bands can be 
considered independently. In this case each band gener- 
ates its own decay rate l/f(g) of the wave function. By 
analogy with Eq. ( p.l6| ) , we can write 



1 1 



In 



E' 



VF2 



(2.19) 
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where is the variance of Uq{p), 



(2.20) 



and Co((?) is the correlator of the averaged potential, 



Co{q) = C{q) [FNiqP)]' 



(2.21) 



Let be the wave vector corresponding to the largest 
^(g), then it is natural to think that — S,{q*)- In other 
words, the localization length should be determined by 
the "optimal band" of harmonics, which we are going to 
find next. 

In view of Eq. ( 2.21 ), two cases have to be distin- 
guished, g* < 2fci? and q* > 2kF- The latter is realized 
for a sufficiently weak "white-noise" random potential, 
the former for the potentials of all other types. 

If g < 2k p, then Co{q, 0) differs fr om C (q, 0) only by a 
pre-exponential factor. Using Eqs. ( 1.12| ) and ( 2.20| ) and 
omitting some unimportant pre-exponential factors, we 
arrive at 



1 1 



In 


E 


{qd)n 




W 


2(3 \ 



(2.22) 



If /3 > 1, then ^(g) given by Eq. (2.22) has the maximum 
at 



d\[3-l 



E 
W 



(2.23) 



Substituting this value into Eq. (2.22) and taking E 
TujJc/2, we obtain Eq. (1.13). 



The qualitative derivation of Eq. ( 1.15 ) goes along the 
same lines. The sole difference is that turns out to 
be close or even larger that 2kF and at the same time 
smaller than 1/d. In this case C'o(g*) is determined by 
Ftq{q^P) rather than by C{q*). (In this case, of course, 
the appropriate width Aq of the bands is much smaller 
than l/d, but the basic idea of dividing the spectrum into 
independent bands stays). 



III. SHORT-RANGE RANDOM POTENTIAL 



In this section we present a more detailed calculation 
of the localization length for the short-range random po- 
tential, d <C As wc m entioned in the preceding section, 
Thoulcss formula ( ^.16| ) is in agreement with the locator 
expansion for Gp, 



GpiO,Mqf;E) 



2t:PE 



nMfc^. ,3,, 



The equivalent of Eq. (3.1) in the general case is 



dq_ 



M 



M=0 ' 
\ n— 1 / m— 1 



2tt 



E 



where pmy stands for 



Pniy = I'^iqi + q2 + ■■■+ Qm-l + T^lm) ■ 



Combining Eqs. ( p.7| ) and (3.2), we arrive at 

oo 



(3.2) 



(3.3) 



(3.4) 



A/=0 



q(M) 



dp 



2t:PE 

n ^0 (9m. P + Pmy) : 



(3.5) 



where S = qJ'^- 

Formula (3.2) is certainly just an approximation. How- 
ever, the model studied in the preceding section showed 
that the the relative error in calculating ^ in this way is 
of the order O {W/E)'^ , which is quite satisfactory. The 
major defect of our approximation is having all the en- 
ergy denominators equal to E. Consequently, this ap- 
proximation scheme does-not capture the phenomenon of 
the resonant tunneling£3 which appears due to anoma- 
lously small energy denominators. Note however that 
our goal is to calculate (In |G(0, y; i?)|). The resonant 
tunneling configurations are exponentially rare and do 
not contribute to this quantity. In the model studied in 
the previous section this can be seen explicitly: the res- 
onant tunneling configurations correspond to E = E' in 
Eq. (2.15) where the integrand diverges. However, the 
divergence is integrable, and moreover has an exponen- 
tially small weight. 

We calculate (In |G(0, y; i?)|) in three steps. First, 
we calculate (|G(0, y; £')p)nr, where the subscript "nr" 
stands for "non-resonant," i.e., with resonant tunneling 
configurations excluded. The reminder of such an ex- 
clusion is essential in this case because even being ex- 
ponentially rare, the resonant tunneling configurations 
yield untypically large |G(0, y; E)\'^ and totally dominate 
the average square modulus (|G(0, y; E)\'^) for sufficiently 
large y (see Ref. |l5| and Ref. |2^) . 

As the next step, we calculate the decay rate of 
{\G{0,y;E)\^)nr defined similarly to Eq. (^, 



- lim ;^ln(|G(0,y;i?)|2)^ 
^2 y^°° 2y ' 



(3.6) 



Finally, we show that £,2 = 

The calculation of (|G(0, y; i?)p)nr can by represented 
with the help of diagrams, one of which is shown in Fig. |^. 
The solid lines in this diagram correspond to the factors 



7 



1 /E, each dashed Une stands for Co (g, p) with appropri- 
ate arguments, and the vortices at the corners bring the 
factors ^n{-)^n{ )- 




FIG. 5. A typical ladder diagram, which describes the 
tunneling in the short-range random potential. 

The diagram shown in Fig. || is of the ladder type. 
It is easy to see that other diagrams (with crossing 
dashed Unes) are negligible. Indeed, consider, for in- 
stance, a diagram where mth and m -I- 1st dashed lines 



of the original ladder diagram are interchanged. This di- 
agram will be proportional to Co{qm, ^Py) where Apj, = 
i,y ^ qnJ^- As discussed in the preceding 
section, the characteristic values of are of the or- 
der of q<t{E), so that the distance Apj, between p„iy and 
Pm+i.y is of the order of q^{E) P. This distance is much 
larger than d, the correlation length of Uq because I ^ d. 
Therefore, Co(gm, ^Py), and thus, the entire diagram are 
small. 

Note also that the neglect of the resonant tunneling 
configurations is assured by omitting the diagrams with 
dashed lines connecting two points of the same solid line 
(either the upper or the lower one). 

The magnitude of the ladder diagram in Fig. || is equal 

to 



dqm 
2ttE^ 



Co (<7m,P2 - Pi) 



(3.7) 



The products of functions $jv in the square brackets 
can be replaced by the integrals over auxiliary variables 
according to the formula 

*»(.+|)*»(.-|)=/^ 



e""' Fr(t,y), (3,8) 



where F]sf{z,y) — Fn[\/ + y^), which enables one to 
obtain a rather simple expression. 



d^pFl^p) Co{pf' e^P'^y/'" 
(27r?2)3 



and finally, 

A very similar calculation gives 

dp 



J7J2M+2 



d'^pF^ip) e'P-y/i'' 



{\GpiO,y;E)\^ 



Jpy/i 



(3.9) 



(3.10) 



(3.11) 



(271/2)2 E^-Co{p)' 

Since the integrands in these formulas oscillate the more 
rapidly the larger y is, the square modulus of the two 
Green's functions de cays with y. Furthermore, compar- 
ing Eqs. ( |3.10 ) and (3.11), we see that the decay rate of 
C and Cp is the same, i.e., that ^2 would not change if 
C were replaced by Gp in Eq. (3.6). In both cases ^2 is 
given by 

6 = I'/v*, (3.12) 
where is the smallest positive root of the equation 

'dq.^dqy 



(2^) 



■Coiq.,qy)e^'''^'" =E\ (3.13) 



Using Eqs. ( 1.12| ) and (2.21), and also the asymptotic 
formula [FAr(^F)]2 ~ l/iirqE^) (vahd for > 1 and 
q kp) we obtain 



dq 



7ri?c(47rgw*)i/2 



g2 

c(o)- 



(3.14) 



If /3 > 1, the integrand has the saddle-point at q = 
[Eq. ( ^.23| )] with the characteristic spread of q around 
being A(j ~ l/{d [ln(i;/M^)](''-2)/2/3}. Using the saddle- 
point method estimate for the integral, and then solving 
the resulting transcendental equation, we obtain 



^InM 
- 1 W 



(3.15) 



As one can see from Eqs. ( ^.12 ) and (3.15), the derivation 
of Eq. (p^ ) will be complete if we demonstrate that ^2 = 
i.e., that ln(|G(0, y; E)\^)^, - (In |C(0, y; E)\^) = o{y). 
(Note that we are interested mainly in the case \E\ — 
fiLDc/l). Since such a calculation is not an easy task, we 
will only show that this relation holds for G^'^\ where 
M K, y/^qj"^). Such M give dominant contribution to 
(|G(0,y;£:)|2)„,., and presumably, to (In |G(0, y; £:)|2) as 
well. 

To average the logarithm, we employ the replica trick, 
(|G|2") - 1 



(ln|G|) = lim 



2n 



(3.16) 



Under the same kind of approximations as above and for 
integer n, (|G(^^)|2")„, is given by 
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(|G( 



M)\2n\ 



n 

m — 1 



2t:PE 



^N{pm)'^N{p 



v) 



M 



nn 

k=l r=l 



.1 (^) 



A- 



(r) 



JO 



(3.17) 



where Em = X^fli 91'"^^^ -Pfc labels the permuta- 



(m), 



tions of the superscripts of the set {g^ ,5^ 
The quantity a'^y stands for P[gf^''^^ + . . . 4 



ig^*" '■'"■']. There are altogether n! permutations Pk for 
each A;; therefore, the complete expression is a rather 
complicated sum of (n!)*^ terms. However, only n\ terms 
in this sum are significant. Indeed, within the adopted 

(r) (r) 

approximation all the terms with — (Jj, \ 3> d for 
at least one of k — 1,2,..., M, should be dropped. It 



'ky 



for the case 
Pm- In this case 
a'j^y = Pr ~ PPi(r) for all k. Therefore, a single 



is easy to find the difference p\.y 
of identical permutations Fi , P2 , • 

Jr) _ Jr) 
Pky 

constraint \pr — Pp^(r) \ ^ d takes care of all the M con- 
straints above. If, on the other hand, some of Pk are not 
the same, then the integration domain is much more re- 
stricted, and the value of the integral is small. Retaining 
only the terms corresponding to identical permutations, 
we immediately find that 

(|GW|2")„,~n!(|G(*'^)|2)J^^.. (3.18) 

Taking the limit 71 0, we obtain from here0 

ln(|G(^^)|2\ -r- 



(ln|G(^-^)|> 



(3.19) 



We think that the same relation holds if G*^^^^ is replaced 
by the total Green's function, G(0, y; E), and so ^2 = 

Finally, let us sketch the derivation of ^ for the "white- 
noise" random potential, d kp^. In this limit one can 
replace Go(q) by C{0)F^{qf) in Eq. (|3T3|) , which leads 
to the following equation on u*. 



cm 

2nP 



F^{2tv,) = E^ 



(3.20) 



The next step is to use the asymptotic formula for 
FN{iy), 

2N+1 




(3.21) 



4i?2, (3.22) 

valid for y 3> fc^^. This way one obtains an approxi- 
mate solution for v^,. Finally, taking \E\ to be hu)c/2, 
one recovers Eq. ( 1.15| ). 



IV. LONG-RANGE RANDOM POTENTIAL: 
DISCRETE LANDAU LEVELS 

In the previous section devoted to short-range random 
potentials, we were able to derive ^ by calculating the 



square modulus of the Green's function. Unfortunately, 
this is not possible for a long-range random potential. 
The physical reason is as follows. 

As we have shown in Sec. the typical distance Apj, 
between the locations of successive scattering events is 
of the order of {P/d)C, where £ is some logarithmic fac- 
tor. If the random potential is long-range, d ^ I, then 
Apy <C d, and so such scattering events can no longer 
be considered uncorrelated. One of consequences is an 
enhanced probability of "pinholes" discussed in Sec. ||. 
In other words, the local decay rate of the wave func- 
tions with distance becomes very nonuniform. In its turn, 
the Green's function G{0,y;E), even with the resonant 
tunneling configurations excluded, exhibits large fluc- 
tuations between different disorder realizations so that 
|ln(|G(0,y;i?)|2)| » |(ln JG(0, ?y; £;)P)n.l, or 6 »e 

One look at Eq. ( 3.17 ) is sufficient to predict that a 
diagrammatic calculation of (In |G(0, y; E)\'^) is bound to 
be very cumbersome. The task is easier within a different 
approximation scheme, the WKB method. Some of the 
formulas corresponding to this approximation have been 
previously wpcked out by Tsukadao and by Mil'nikov 
and Sokolov.y 

Suppose we want to find the solution of the 
Schroedinger equation Uo(f>{py) — E4>{py). Let us seek 
the solution in the form 4>{Pv) = exp[?S'(/3j,)]. The action 
S{py) can be expanded in the series of the small parame- 
ter l/d. In the lowest approximation, S{py) must satisfy 
the Hamilton- Jacoby equation. 



,2 dS 



so that 



HPy) ^ exp 



Py 



j2 I dvPxiv) 



where Px{Py) is a solution (generally speaking, a complex 
one) of the equation 



Ua{p^,Py) = E. 



(4.1) 



The meaning of this equation is quite transparent. It is 
known that the motion of the guiding center in classically 
permitted regions is a drift along the level line s of the av- 
eraged potential Uq (see Ref. ^l|). Equation (4T) means 
that the trajectory of the guiding center in classically 
forbidden regions is still a level line although obtained 
by analytical continuation to the three-dimensional space 
{(w, V, Py)} where u = Re px and v = — Im p^- 

The WKB-type formula for the Green's function in the 
guiding center representation is. 
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Gp{pi,p2\E) 



sgn(p2 - Pi) 



V L J o„=ni L 



exp 



dpy p'x'Hpv) 



(4.2) 



Py=P2 



(n) 

where the super scri pt n labels different solutions px = 
Un — ivn of Eq. (O) in the complex half-space sgnu„ = 
sgn(p2 — Pi)- If we study the tunneling from point (0, 0) 
to (0, y) where y > 0, then this will t ypic ally b e th e up- 
per half-space, w„ > 0. Using Eqs. ( ^.7[ ) and ( |4.2| ) and 
neglecting all the pre-exponential factors, we obtain the 
following estimate for G(0, y; £'), 



exp 



UnW+ul{y) I 



2/2 



p / dpy P^'iPy) 



Since Uq is a random potential, it is natural to assume 
that the phase factors corresponding to different n in this 
sum are uncorrelated; therefore, 



exp 



pi dpy Vn{Py) 



Consequently, ^ can be calculated as follows. 



1 ■ r 1 
— — mm lim — 



uliy) 



+ / dpy VniPy) 



(4.3) 



(4.4) 



It is possible to demonstrate that the first term in the 
square brackets is typically much smaller than the sec- 
ond one, which leads to 



/2 



(4.5) 



where (vn) is the average "height" of the nth level line. 



{vn) = lim - / dpy VniPy). 
y^oo y ' 



The problem of calculating min (u„) turns to be rather 

n 

difficult and we have not been able to solve it exactly. 
However, we will give arguments that min (w„) ~ [w, 

n 

was introduced in the previous sectio n, se e Eq. (3.13)]. 
Therefore, ^ is still given by formula ( 1.13 ). A more ac- 
curate statement is as follows. If the individual Landau 
levels are well resolved in the density of states, then for 
arbitrary range d of the random potential f can be found 
from the same "master" equation 



(4.6) 



Let us familiarize ourselves with the properties of the 
equipotential contours (level lines) of the potential Uq in 
the half-space {{u, v, Py),v > 0}. An important property 
is the contour density P{v, E), 

P{v, E) = '^(« - '^(" ~ "»)^ 

(for isotropic random potential Uq this quantity does not 
depend neither on u nor on py). Function P{v, E) proves 
to be the sum of three terms. 



P[v, E) = Pi[v, E) + P2{v, E) + P3(S) 5{v), 
where 

E^S^ 



Piiv,E) = 
P2[v,E) 



1 



2ny/QR 

and Q, R, S, and T are as follows, 

Q = Co(2w,0) + Co(0,0), 

i? = Co(2zt;,0)-Co(0,0), 

S=\dQ/dv, T^\d^Q/dv 



(4.7) 

(4.8) 
(4.9) 



(4.10) 



Except for the term P3(v) our formulas are in agreement 
with Refs. |l^ and This term, however, does not play 
any role in the subsequent calculation. The derivation of 
Eqs. (g^ and (g||) can be found in the Appendix. 




V (d) 
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FIG. 6. Functions Pi{v,E) and P2{v,E). Vertical axis 
is in units of . The parameters used in generating the 
plot are /3 = 2 and E/W = 5. 

The behavior of functions Pi and P2 is illustrated by 
Fig. |6| Function Pi has a sharp maximum at u = 
where Q ~ E^. Away from the maximum it is expo- 
nentially small. Function P2 is exponentially small at 
f < w*, and assumes the asymptotic form P2(v,E) cx 
^(2-/3)/(/3-i) g^i- ^ ^ Yor example, if /3 = 2, then 
P2{v,E) l/in(f) (see Fig. |). The ratio P1/P2 eval- 
uated at the point is of the order of In ^ 1 in 
the case of interest. At v close to our expression essen- 
tially coincides with Eq. (5.20) of Ref. Note however, 
that the latter equation is off by An. 

Let us clarify the origin of the sharp maximum in 
P{v,E) a,t V = v^,. To this end the concept of bands of 
harmonics introduced in Sec. || is very helpful. So, let us 
consider a band with in the interval (g— 5 Ag, q + ^Aq) 
where q ^ Aq > 0. The amplitude of each harmonic 
Uo{qx,Py) is enhanced by the factor e''"'" upon the ana- 
lytic continuation into the upper half-space. Therefore, 
the typical val ue of t he combined amplitude Uq{p) of the 
band [see Eq. (2.18) for definition] grows from Wq{Q) ^ 
[Co(g,0)Ag]i/2 at « = to Wq{v) ^ [Co(g, 0)A9]i/2 e?" 
sX V > 0. Being the product of the rapidly decreasing 
function Co{q,0) and the exponentially growing factor 
e'", this quantity has a sharp maximum at q = q*iv), 



(4.11) 



provided that /3 > 1 and v ^ d. [In view of Eq. ( p.l5[ ), 
this formula is just another parametrization of origi- 
nally defined by Eq. (2.23) as a function of E\. Gener- 
ally speaking, the width of the maximum depends on v 
and /? but in a particular example, /3 = 2, it is simply 
Aq ^ 2- Hence, at a given "height" v in our three- 
dimensional space, the potential Uo{px,Py) is typically 
dominated by the band of harmonics of width Aq cen- 
tered at q*{v). Consequently, the spatial dependence of 
Uq is almost plain-wave-like, Uo{px,Py) (x e.'^P^i*'^^) , This 
prompts the decomposition 



Uo{p) - 



V{p) 



(4.12) 



of C/q into the "oscillating part" e^p^i'('") and a "smooth 
part" V{p). 

The intersection points of the level lines Uq = E with 
a vertical plane Py = i] satisfy the system of equations 



-^InM, 

g,(w„) E ' 

1 , , 27rn 



arg V 



(4.13) 
(4.14) 



The modulus | V| of the complex Gaussian variable V has 
the Maxwellian distribution with the characteristic width 
W given by 



W'^{\Vip)\' 



-2vq.M [v) 



Q, 



(4.15) 



where Q is defined by Eq. (4.10). In the first approxi- 
mation let us neglect the dependence of V on p^, (which 
corresponds to the limit of infinitesimally narrow band, 
Aq 0), then P{v, E) is given by 



P{v,E) 



q* 
2^ 



md\v\ 

W2 



V q* 



In 



Doing the integration, we arrive at 
P(.,i.).^e-VQ, 



(4.16) 



which, in fact, coincides with Pi(v,E) provided that 
V ^ d. In this approximation Un form an equidistant 
set, Un+i — Un — 27r/(7* <C 1/d, while Vn does not depend 
on n. The equipotential contours resemble a number of 
uniformly spaced parallel rods, which "soar" above the 
real plane v — 0, staying very close to the "standard 
height" u = for most values of py. Indeed, as dis- 
cussed above (see also Fig. H), function Pi has a sharp 
maximum at w = . Since c? In Q /dv ~ 2q^ , the width of 
the maximum is of the order of 1/ q*. This is, of course, 
can be seen directly from Eq. ( 4.13 ): if |V| has its typical 
value, W, then ti„ = w*. Fluctuations of [V"! change the 
logarithm (typically) by a number of the order of unity; 
therefore, ordinarily |w„ — w^l ~ 1/9*. Significant devi- 
ations from the "standard height" are exponentially 
rare. 

We believe that such a description accurately portrays 
the behavior of the relevant equipotential contours in 
the upper half-space, which means that ~ for 
such c onto urs, and therefore, 1^ is given by the old for- 
mula ( 1.13| ). 

As one can see from Eqs. ( p77[ - [1.9[ ) and (4.16), our ap- 
proximate treatment does not capture the term P2{v, E) 
in P{v,E), which seems to be important for v > w*. 
Therefore, the possibility of min„(w„) being larger than 

can not be totally ignored. Although we can not rig- 
orously prove that min„(t;„) = w*, we managed to find 
the upper bound for min„(i;„). 



l{Vn) < W* + 0{l)d, 



(4.17) 



based on the following percolation- type arguments. 
Suppose that there exists a level line number to, 
,Py] = E, which is totally contained in a slab 
{{u,v, py),0 <v< vq}. In this case min„(i;„) < (vm) < 
vq. This prompts considering the following percolation 
problem. Let us call "wet" all the points {u,v,py) of the 
slab, which satisfy the conditions 

Re Uq{u — iv, Py) < E, Im Uq{u — iv, py) < 0. 
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As the thickness vq of the slab increases from zero, it 
should eventually reach a critical value Vc, at which the 
percolation through "wet" regions first appears. The per- 
colation threshold Vc is at the same time the upper bound 
for min(?;„). 

n 

Obviously, no percolation exists for vq < when 
the wet regions occupy exponentially small fraction of 
the volume. On the other hand, if vq S> w*, then 
((ReC/o)^) - ((Imf/o)^) - Q/2 > E^, and approxi- 
mately a quarter of the entire volume is wet. This greatly 
exceeds the volume fraction 0.17 required for the-perco- 
lation in continual three-dimensional problems;t3 thus, 
such Vq are high above the percolation threshold. Fur- 
thermore, the relation betweeiilhe percolation thresholds 
of the continuum and of a filrrO suggests that m ust b e 
equal to v^, +Cd, where C ^ 1, which leads to Eq. (4.17). 

Finally, let us comment on the relation of our approach 
to that of Raikh and Shahbazyanllj already mentioned 
above. In Ref. ^ the idea of the complex trajectories 

Px^\py) — Un{py) — ivn{py), which is the basis of our cal- 
culation of the Green's function, was introduced. How- 
ever, an approximation early in their analysis [dropping 
of the cross-product in Eq. (2.8) of Ref. |lj] led the au- 
thors of Ref. |lj to the result, which in our notations can 
be written as follows: 



Gp(0, y;E) = ^Y. [5C/o/ap(")]-i exp 



P 



\G\' 



E '^^P { - ^ K {y/2) + yv4y/2)] 



instead of Eq. (4.3) and 



ay) " 



\{ym 

y 



■Vniy/2) 



(4.18) 



instead of Eq. (|^. If y < d exp{E^ /W^), the mini- 
mum is typically supplied by one of the trajectories whose 
"height" Vn is not too different from the standard value 
of u*, |w„ — w*| ^ d. This range of w„ allows for m„ in 
the range |u„| < \/yd. Therefore, the optimal trajec- 
tory is typically one of about M = \fyjd trajectories 
closest to the x = axis. Note that M increases with 
y. Eventually, at ?/ ~ d exp(iJ^/WQ ) it becomes expo- 
nentially large so that there is an appreciable probability 
that one out of such M trajectories has Vniy/'^) < w*/2 
and |u„(?//2)p < yd/2. The decay rate l/^(y) in this case 
is determined by such an untypical traject ory, a nd so ^(y) 
is significantly larger than our result, Eq. (|l1^ ). Further- 
more, at y 3> d exp(£'^/T/FQ ), there is a finite probability 
of finding u„ { y/2) exactly equal to zero [due to the third 
term in Eq. (4.7)]. Therefore, in the asymptotijc limit 
y ^ CO the method of Raikh and ShahbazyanO yields 
a very surprising result 1/^(2/) — *■ 0. As we explained 
above, this is a consequence of an effective substitution 
of the original two-dimensional random potential by a 
potential depending on a single coordinate. 



(the derivative is taken at py = y/2). As one can see, 
in their method Gp(0, y; E) is determined not by t he en- 
tire tunneling trajectory, as it should [see Eq. (4.2)], but 
only by its midpoint py — y/2. As a result, the point- 
to-point Green's function G(0, y; E) obtained by their 
method starts to be in error strictly speaking already 
for y ^ d. ThCjjmly case where the method of Raikh 
and Shahbazyantil works for arbitrary y is the case of 
a one-dimensional potential, e.g., Uo{x,y) = Uo{x). A 
potential of this type has the same magnitude at the 
midpoint y/2 and at all other points y. (From the formal 
point of view, only for this type of potentials the afore- 
mentioned cross-product vanishes for any y and therefore 
can be dropped). 

The reduction of the properties of the entire trajec- 
tory to the properties of a single point destroys the self- 
averaging of the localization length, which is-a well estab- 
lished property of other disordered systemsB3 Because of 
that, it becomes necessary to treat ^ as a function of y. 
As we mentioned above, for a given disorder realization 
the value of G(0, y; E) given by the formulas of Ref. |lj 
is in error already for y ^ d. However, this is not so for 
the quantity (In |G(0, y; -E)|). The logarithmic averaging 
manages to mask the error so that ^(y) re mains close to 
the correct asymptotic value given by Eq. ( jTrsI) for suf 
ficiently short distances, y <C d ex.p{E'^ /Wq). Indeed ^ 
one uses the last equation above, then one obtains 



V. OVERLAPPING LANDAU LEVELS 



if 



This section is devoted to the derivation of Eq. ( 1.10 ). 
An important difference from all the preceding calcula- 
tions is that the potential energy E = Ep — 'hujc{N+^) is 
smaller than the amplitude Wq of the averaged potential 
Uo{p). In this case the level lines J7o[pi"\ Py] = E intro- 
duced in Sec. [V, stay predominantly in the real plan e 
V — 0. [Their density is given by P3{v,E), see Eq. (4.7)]. 
However, for E ^ Q the percolation of the level lines in 
the y-direction discussed in the end of the previous sec- 
tion, still requires brief excursions into the upper complex 
half-space. Such excursion link together the large closed 
loop contours = E oi typical diameter ^porc(£')- The 
localization length ^ is still given by Eq. ( |4.5| ), which 
leads to the estimate 



l"^ Cpcrc(-E) 



Ap 



(5.1) 



where the subscript "ex" indicates that the averaging is 
performed only over the "excursions" , i.e., over the parts 
with u > 0, with Ap being the typical length of such 
parts. Similar to the caseJV = discussed previously 
by Mil'nikov and Sokolovjl3 (z;)ox ~ Ap - {\E\/Ulif/'^ 
where C/q ~ W^/d^ is the typical value of the second 
derivatives of the averaged potential Uq; therefore, 
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Ed? 



(5.2) 



The calculation of tho-quantity ^porc(£') is a subject of 
statistical topography,ll3 and the following results are es- 
tablished. Denote by C\ the integral, 



Ca = 



i<|<?|A<l 



d'^q 
(2^ 



Coiq). 



For slowly decaying correlators CQ{q), such that C\ ~ 
X^^^ with H < 3/4, ^porc is given by (besides the re- 
view, Ref. ^ see the original works, Ref. ^ ), 

W(£;) ^dlW^o/^r/''. (5.3) 
Otherwise, i.e., if Cx decays faster than A~'^/^, thenEl 

^perc{E)^d\Wo/E\^/\ (5.4) 




W/E 



FIG. 7. The percolation length ^porc as a function of 
the energy E (schematically). The labels "2" and "4/3" 
indicate the power-law exponents in the corresponding in- 
tervals, see Eq. ([3.5|). 

Let us now determine the conditions under which for- 
mulas (pT^) and (5.4) become applicable. The issue is 



complicated by the fact that Co{q) is the product of two 
terms, Co{q) and [FN{ql^)f, see Eq. (|3l]). The for- 
mer decays exponentially starting from q ^ 1/d. The 
latter remains close to one at q ^ ^/Rc, then be- 
haves according to a power-law, [Fjv((jP)]^ ~ 1/q for 
1/Rc < q < 2kp, and finally, decays exponentially at 
q > 2kF- Such a complicated behavior results into three 
different regimes (l.lOa-c). 



The simplest is the case Re < d, where \FN{qr)V 1 



for all relevant q. In this case Eq. (5.4) applies and 
also Wq ~ W. The localization length is given by 



Eq. (1.10a), which coincides with the result of Mil'nikov 
and Sokolov.y 

If Rc S> d, then the situation is more complicated. 
In this case Cx is proportional to 1/A, i.e., H ~ ^ ioi 
d <C A <C Rc, yet decays faster than for X > Rc- 

As a result, both Eq. (5.3) and Eq. (5.4) for ^pcrc may 
apply, depending on E, 



£,pcrc{E) ^ < 



d^ 

Rc 
Rc 



W 

'e 

d W 

Tc^ 



d \E\ 
— < - — - < 
Rc W 



d 

Rc ' 



(5.5) 



\E\<W—, 



which is illustrated by Fig. W. Combining Eq. (5.2) where 

and using \E\ = fiLOcjl, 



Wq - W^fdlRc with Eq, 



one obtains Eqs. (1.10b) and 



z). 



VI. DISCUSSION OF A MORE REALISTIC 
MODEL AND ITS COMPARISON WITH THE 
EXPERIMENT 

To make the connection with the experimental practice 
we will consider the model where the random potential 
is created by randomly positioned ionized donors with 
two-dimensional density rii set back from the electron 
gas by an undoped layer of width d. We will assume that 
1/d^ <C Tii <C v?(P . In zero magnetic field the random 
potential can be considered a weak Gaussian random po- 
tential with the correlator 



C(g) = Tp'ni[<P'aB) 



2 g-29d 



(6.1) 



where as = fi^K/me^ is the effective Bohr radius (see Ap- 
pendix B of Ref. 11). This formula corresponds to /3 = 1 
and d replaced by 2c? in Eq. (1.12). Deriving Eq. ( |6.lD , we 
took into account the screening of the donors' potential 
by the electron gas described by the dielectric functions 



K 1 



asq 



q < 2kF. 



(6.2) 



This model remains to be accurate in sufficiently weak 
fields, where the Landau levels overlap and the density 
of states is almost uniform, like in zero field. In stronger 
fields (the boot-shaped region AHIGEB in Fig. g), the 
density of states develops sharp peaks at the Landau 
level centers separated by deep minima. This strongly 
influences the property of the electron gas to screen the 
external impurity potential. Different aspects of such a 
scree ning in weak magnetic fields have been addressed in 
Refs. ^9|-^. The screening can be be both linear and 
nonlinear, depending on the wave vector q. 

The concept of nonlinear screening has been developed 
by SUclovskii and Efros iaitially for the three-dimensional 
case.El Gergel and Suris£3 have extended it to the two- 
dimensional case. The influence of a strong magnetic field 
on the nonlinear screening has been studied in Ref. ^ 
and especially in Ref. 

Nonlinear screening is realized for sufficiently small q, 
q < ^noni, and is enhanced compared to the zero field 
case [Eq. ( |6.2[ )]. The threshold wave vectpi; (?noni is a 
complicated function of the magnetic field.Eil If 9noni is 
smaller than 1/Rc, there exists an intermediate range of 
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Q, 9noni < 9 < 1/-Rc, where the screening remains i n th e 
Hnear regime but is suppressed compared to Eq,-|d6f2|) 
The corresponding dielectric function is given bycSE 



eiq) ^ ^ 1 + ^ 



E^q 



(6.3) 



At even larger q, q > max{(7noni, there is no es- 

sential change in the screening properties brought about 
by the magnetic field. 

Clearly, in such magnetic fields the model of a Gaussian 
random potential with the correlator (xl) becomes an 
oversimplification. Fortunately, the localization length 
^ should not be strongly affected by this. Indeed, ^ is 
sensitive only to the combined amplitude of the narrow 
band of harmonics with wave vectors g ~ q*. It can be 
shown that for /3 equal to one, ~ 2kF (just like for 
the "white-noise" potential). This wave vector belongs 
to the last group of q for which there is no change in the 
screening properties. In fact, at such q the screening is 
equally ineffective both in zero and in arbitrary strong 
magnetic field, e(2fcF) — k. 

The parameter kpd is almost always larger than one in 
the experiment; t heref ore, in strong fields is expected to 
be given by Eq. ( 1.13 ). For (3—1 this formula becomes 
remarkably simple. 



d 



Rc 
kpd 



The simplest way to obtain this expression is to start 
with Eq. ( 1.13 ) for /3 > 1, take the limit /? ^ 1, and then 
make a replacement d ^ 2d (see above). For kpd close 
to unity, however, one should use a refined formula 



Rc 



kpd - 



(6.4) 



which follo ws from the master equation ( |4.6| ). For- 
mula ( |6.4| ) should hold in the boot-shaped region 
AHIGEB of Fig. |l|. The corresponding range of N can be 
expressed in terms of the dimensionless parameters kpd 
and rii/n. We will concentrate on the ca se n j ~ n, which 
we call "standard." In this case Eq. (6.4) is valid for 
N 5, 2k pd. A t larger N, the localization length is given 



by Eq. ( |1.10[ ) so that the dependence of ^ on is super- 
linear. Note that the "standard" case corresponds to the 
straight-line "trajectory" passing through the points E 
and F and s hown by the arrows in Fig. |l|. Between these 
points, Eq. (1.10) reduces to 



Tc ^ {kpdfi^ ' 



2kpd <N< C{kpd) 



5/3 



(6.5) 



where C is some undetermined numerical factor. At even 
larger N, ^ becomes exponentially large, which precludes 
its accurate measurement at experimentally accessible 
temperatures. For this reason we do not give the explicit 
formula for the localization length at such N. 



Unfortunately, the published experimental data on the 
low-temperature magnetoresistance away from the QHE 
peaks is limited to the measurements done by Ebert et 
al.U more than a decade ago. In most of their samples d 
was equal to 6 nm, which corresponds to kpd — 0.9. On 
the one hand, this places the samples close to the line 
HI in Fig. [l[ On the other hand, the "standard" case 
corresponds to the line EE. In fact, there is no contradic- 
tion here because in the "standard" case Ep/W ~ kpd. 
Therefore, the as kpd tends to unity, the boot-shaped 
region AHIGEB shrinks and the lines EF and HI become 
quite close. 

With the help of Eq. (L4), we converted Tg reported 
by Ebert et alQ into ^. Only even ly w ere s elected for this 
analysis. The constant factor in Eq. ( |1.4D was chosen to 
be 6.0 following Nguyen.t3 For filling factors v = 2 and 
ly = A, i.e., N = 1 and N ~ 2, there is a good agreement 



with Eq. (6.4) if the logarithm in the denominator is ne- 
glected. At larger J^, 6 < < 12, empirical ^ behaves 
roughly like 



^ « 2.2N^^^Rc 



(6.6) 



Formally, this is in agreement with Eq. (6.5). However, 
since the parameter kpd is close to one, the random po- 
tential can not be treated as long-range, and so the per- 
colation picture used in derivation of Eq. (6.5) is not 
quit e ju stified. Perhaps, it is more resonable to treat 
Eq. ( |6.6[ ) as an approximation to the exponential depen- 
dence, Eq. (L8), within a limited range of ly. 

Another comment is in order here. Formula ( |l.4| ) is 
based on the assumption that the interaction energy of 
two quasiparticles separated by a typical hopping dis- 
tance r is given by the Coulomb law. As was s how n 
in Ref. with a dielectric function given by Eq. (|6.3|), 
the Coulomb law is realized only at sufficiently large dis- 
tances, r ^ R^/qb- Large hopping distances r corre- 
spond to low temperatures, T = (^/4r)^To. Therefore, 
formula ( |l.4| ) is expected to hold only for 



(6.7) 



Experimentally, such low temperatures may be hard to 
reach, especially in weak fields where Rc is large. There- 
fore, we think that a different type of measurement may 
be more promising. One has to measure both the tem- 
perature and the current dependence of the magnetore- 
sistance. Comparing the results of these two types of 
measurements, one can determine the effective temper- 
ature Toff for each value of the current density j. As 
discussed in Ref. ^, the relation between Tcs and j is as 
follows, 



ks Tcffij) ~ 0.5 ejpxy 



(6.8) 



where pxy is the Hall resistivity. Equation (6.8) does not 



involve the dielectric function and therefore, is expected 
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to hold even when Eq. (1.3) and the corresponding cur- 
rent dependencJ^ p^x oc e""^/"^''/'^"''''^^-' do not match the 
magnetoresistance data exactly. 
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APPENDIX A: LEVEL LINE DENSITY P{v,E) 

The equation Uq = E is equivalent to the system of two 
equations R = E and / = 0, where R = Re Uq{u — iv, py) 
and / — — ImJ7o(u — iv,py). To calculate P{v,E) we 
need to know the Jacobian J = \d{R, I)/d{u,v)\. With 
the help of the Cauchy relations, 

dR/dv = dl/du, dR/du = -dl/dv, 

the Jacobian can be written as 



J Ry + ly , 



where subscripts denote the partial derivatives. There- 
fore, for all V ^ 0, 



As expected, PsiE) is exponentially small for \E\ ^ Wq 



X exp ( — -v^K 



(Al) 



where v is the four-component vector (i?, Ry , /, ly )^ and 
K is its autocorrelation matrix, K^n = {vmVn)- Matrix 
K turns out to be block-diagonal, the blocks being 2x2 
mat rices . Matrix elements of K can be calculated using 
Eq. (1.5). For example, the upper diagonal block has ele- 
mcntsKii = Q/2, K12 = K21 = S/2, and K22 = R/2 [see 
Eq. (4.10) for definitions]. The subsequent integration 
over Ry and ly in formula ( |Al[) becomes a trivial task 



and yields the result represented by Eqs. (4.8) and (4.9). 

The case v = requires a special consideration because 
the condition / = is satisfied identically. It is easy to 
see that a level line would typically branch upon the in- 
tersection with the V = plane. The additional branch or 
branches are totally contained in this plane, which gives 
the delta-function-like contribution to function P{v,E). 
The coefficient in front of the delta-function is given by 
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